METHOD FOR REDUCING NOISE, AND 
PULSE PHOTOMETER USING THE METHOD 



BACKGROUND OF THE INVENTION 

The invention relates to a devioe for measuring physiological 
information through use of photoplethysmograms measured at a plurality of 
wavelengths, and. more particularly, to an improvement in .a pulse photometer 
used for diagnosing a problem in a respiratory system or a drcutatory system. 

A known technique for measuring a concentration of light absorbing 
materials contained in blood during medical care is pulse photometry intended 
for measuring an oxygen saturation SpOa, the concentration of abnormal 
hemoglobin sudi as cart»oxyhemogtobin (COHb) and methemogiobin (MetHb), 
or a concentration of injected pigment A device for measuring oxygen 
saturation SpO^ in particular is called a pulse oximeter. 

The principle of the pulse oximeter is to irradiate onto living tissue 
light rays of a plurality of wavelengths exhibiting different light absorbing 
characteristics with respect to a material of interest and to determine the 
concentration of the material on the basis of a pulse wave produced through 
consecutive measurement of quantity of transmitted Jight. 

Japanese Patent No. 3270917 (of., Claims 1, 2 and Figs. 2 and 4) 
disdoses detemiination of a saturation level of oxygen in arterial blood and a 
concentration of light-absorbing materials in the same by: emitting two 
light-rays of different wavelengths onto living tissue; plotting graphs with levels 
of the resultant two pulse waves originating from transmitted light being taken 
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as vertical and horizontal axes, respectiv ly, to thereby determine regression 
lines; and determining the oxygen saturation and the concentration of 
light-absorbing materials in the arterial blood on the toasts of gradients of the 
regression lines. This technique enhances measurement accuracy and 
reduces power consumption. 

However, a large amount of computing operation is still required to 
determine regression lines and gradients thereof through use of many sampled 
data sets pertaining to pulse waves of respective wavelengths. 

The principle will now be descrit>ed by taking, as an example, a pulse 
oximeter for measuring the oxygen saturation In arterial blood. 

This technique is not limited to a pulse oximeter but Is also applicable 
to a device (or a putse photometer) which, on the basis of the principle of the 
pulse photometry, measures abnormal hemoglobin (carboxyhemoglobin or 
methemoglobin) or a light absorbing material in blood such as pigment Injected 
in blood. 

Fig. 1 is a Not* diagram of such a pulse oximeter. The device will 
now be described by taking the pulse oximeter as an example. 

At the time of measurement, a probe 1 is usually attached to an area 
subjected to the measurement. 

A processor 8 produces light emission timings for a light emitter 
(LEDs) 2 which procfcjce light rays of different wavelengths, that is, a red-light 
emitting diode (hereinafter abbreviated "R-LECT) and an infrared-light-ray 
Gmltting diode (hereinafter abbreviated NR-LED"). The R-LED and the 
IR-LED, both being provided in the probe 1 , are driven by an LED driver 4, to 
therelsy alternately emit light. 



The (ight emitt d by the LEDs 2 passes through the area subjected to 
the measurement to which the probe is attached, and is received by a light 
receiver (photo diode; PD) 3 provided in the probe 1 . 

The thus-TOcelved signal is converted from light into electricity by the 
light receiver 3 and is subsequently converted into a voltage by an input 
section 5. 

Components reflecting optical characteristics of pulsation of the area 
subjected , to the measurement appear in this signal in the form of AC 
components. 

By a demodulator 6, the signal output from the light receiver (PD) 3 is 
separated and demodulated Into a waveform of an infrared light ray (IR) and 
that of a red fight ray (R). The waveforms are digitized into signal trains by an 
A/D converter 7. The signal trains are sent to the processor 8, where the 
signal trains are processed, thereby computing, e.g., oxygen saturation Sp02, 

The signal trains that have been digitized by the AID converter 7 
corresponding to infrared light (IR) and red light (R) form respective measured 
data sets. 

Light absorbance A of the light that has passed through tissue is 
expressed by Equation (1), with the attenuation of light in Wood being taken as 
Ab and the attenuation of light in other tissue being taken as At 

A^Ab+At^ Jniliri) - ]n(Ioui) = ECD +At (1 ) 

wherein: 

lout: Intensity of transmitted light 



iin: Intensity of incident light 

E: Absorption coefficient of blood 

C: Concentration of hemoglobin In blood 

O: Thickness of blood 
When the blood thickness has assumed D+AD as a result of pulsation of blood, 
Equation 1 is transformed as follows. 

A-i-Mb-bi(Iw)-ln(Iout)+AhCImi) = £C(D^ .(2) 
Subtracting Equation 1 from Equation 2, Equation 3 is obtained. 

AAb = -A ln(/oi//) = ECM) (3) 

Equation 3 is defined for wavelengths A,1 , X2 as follows. 

AAt\ = -Aln(/cwi/,) ^ £;CAD (4) 
AAbj^ = -Alii(/cw/j) =: E^CAD (5) 

Subscripts 1, 2 indicate that the con^espondlng terms are relevant to 
wavelengths Xi, Xi. Provided that Equation 4/Equatlon 5 = <&s, we have 



a> is determined by measuring Aln(loutl), Aln(lout2), whereby the 
oxygen saturation Is determined as follows. 



Sp02 = f(<I>fi) 



(7) 



When the Influence of body movement or a like factor other than the 
component reflecting pulsation is superimposed on Aln(lout1}, Aln((out2), tight 
absorbance of each of the wavelengths is defined as follows. 

A4«l =~Alii(Aw/;)=^CAD+£>?iC«ADw (8) 
A4*;=i-Aln(/o«/;) = £,CAD+£>i,OiAD/i (9) . 

wherein: 

En: Absorption coefficient of the noise source 

Cn: Concentration of the noise source 

Dn: Change in thicl<nes5 of the noise source 
A waveform (pulse wave) reflecting a pulsation component is not solely 
obtained, and a signal on which noise is superimposed is observed. In this 
case, a light absort>ance ratio 4»' is defined as follows: 

Mt\ ^ li\n(Jaui\) E,CAD+En,CnADn 

A4*; Aln(7cw4)"'£jCAZ>+£>iiCwA2>i ^ ^ 

Hence, the light absort>ance ratio does not coincide with the oxygen saturation 
in arterial blood. 

Fig. 2A shows waveforms which have been subjected to: processing 
in which transmitted light intensity data at respective wavelengths at 
predetermined time intervals are logarithmic-computed; and then a mean of 
the logarithmic>-computed values is made to zero or high-pass filtering is 



performed with respect to the logarithmic-computed valu s. 

Fig. 2B shows a graph through use of waveform data measured 
substantially at the same time when the waveforms shown In Fig. 2A are 
measured (the time is sufficiently shorter than the duration of a frequency 
component of a pulse wave) with the amplitude of infrared light being taken as 
a horizontal axis and the amplitude of red light being taken as a vertical axis. 

If the observed data reflect a mere pulse wave component, the graph 
essentially assumes the shape of a straight line. The gradient of the straight 
line shows the light absorbance ratio a»s. 

However, when noise is superimposed on the data as shown in Fig. 
3A, there Is obtained a result Into which the light absorbance ratio of noise and 
the light absorbance ratio of pulsation are merged as shown in Fig. 3B, in 
contrast with the case of a pulse wave component. 

As mentioned above, when measured pulse wave data include noise, 
an accurate light absorption ratio cannot be measured, and hence noise must 
be reduced. 

Hitherto-known techniques for redudng this noise are a frequency 
analysis technique and an independent component analysis technique. 

However, when the frequency of noise overlaps the fundamental 
wave of a signal component of measured pulse wave data or a harmonic wave 
of the same, the frequency analysis technique encounters difficulty in 
extracting a signal component. 

According to the independent component analysis technique, an 
observed signal cannot be dissolved Into Independent components beyond 
measured waveform data. When a plurality of noise sources are present. 



6 



difficulty is encountered in dissolving the data. 

SUMMARY OF THE INVENTION 

It Is therefore an object of the invention to provide a method for 
determining a light absorbance ratio by readily reducing noise from an 
observation signal on which noise is superimposed, to thereby extract pulse 
wave data. . ... 

It is also an object of the invention to provide a pulse photometer 
using pulse wave data from which noise has been reduced by the above 
method. 

In order to achieve the above objects, according to the invention, 
there is provided a method of processing an observed pulse wave data, 
comprising steps of: 

irradiating a Tmng body with a first light beam having a first 
wavelength and a second light beam having a second \A^velength which is 
different from the first wavelength; 

converting the first light beam and the second light beam, which have 
been reflected or transmitted from the living body, into a first electric signal 
corresponding to the first wavelengtii and a second electric signal 
corresponding to the second wavelength, as the observed pulse data; 

computing a light absorbance ratio obtained from the first electric 
signal and the second electric signal, for each one of frequency ranges dividing 
an observed frequency band\ and 

determining that noise is not mixed into the observed pulse wave data 



in a case where a substantial match exists among light absorbance ratios 
comput d forth r spectiva frequ ncy ranges. 

Preferably, the xlstenc of th_ substantial match of the light 
absorbance ratios is determined with regard to frequency ranges in which at 
least one of the first electric signal and the second electric signal has relatively 
large powers. 

According to the invention, there is also provided a method of 
processing an observed pulse wave data, comprising steps.of: 

irradiating a living body with a first light beam having a first 
wavelength and a second light beam having a second wavelength which is 
different from the first wavelength; 

converting the first light beam and the second light beam, which have 
been reflected or transmitted from the living body, into a first electric signal 
conresponding to the first wavelength and a second electric signal 
corresponding to the second wavelength, as the observed pulse data; and 

whitening the first electric signal and the second electric signal by an 
affine transformation using a known light absorbance ratio, in order to separate 
a pulse signal component and a noise component wtiich are contained in the 
observed pulse data. 

Preferat>iy, the affine transformation Is performed with the following 
equation: 

Where, S Is the pulse signal component, N is the noise component, s1 Is the 
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first electric signal, s2 is the second electric signal, ^ = Tan"^<&, * is the known 
light absorbance ratio, and B is a value selected from a range of ^ to (n/2 - 4>), 
and 

wherein 6 is so selected as to make a norm of the pulse signal 
component minimum. ^ 

Here, it is preferable that the signal processing method further 
comprises steps of; 

computing a light absori^ance ratio obtained from the first, electric 
signal and the second electric signal, for each one of frequency ranges dividing 
an observed frequency band; and 

determining that noise is not mixed into the ot>served pulse wave data 
in a case where a substantial match exists among light absorbance ratios 
computed for the respective frequency ranges. 

Here, one of the light absorbance ratios, which are determined that 
the noise is not mixed therein, is used as the kriown light absorbance ratio. 

Preferably, the signal processing method further comprises a step of 
obtaining a slgnal-to-nojse ratio of the observed pulse wave data by performing 
a frequency analysis with respect to the pulse signal compcHient and the noise 
component at each of predetermined frequencies. 

Preferably,, the signal processing method further comprises a step of 
displaying a pulse wave of the living t)ody, based on the pulse signal 
component 

Preferably, the signal processing method furtiier comprises a step of 
calculating a pulse rate of the living body based on the pulse signal 
component 
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According to the invention, there is also provided a method of 
processing an observ d pulse wave data, comprising steps of: 

irradiating a living body with a first light beam having a first 
wavelength and a second light beam having a second wavelength which is 
different from the first wavelength; 

converting the first light beam and the second light beam, which have 
been reflected or transmitted from the living body, into a first elecUic signal 
corresponding to the first, wavelength and a second .electric signal 
corresponding to the second wavelength, as the observed pulse data; 

whitening the first electric signal and the second electric signal to 
separate a pulse signal component and a noise component which are 
contained in the observed pulse data, for each one of frequency ranges 
dividing an observed frequency band. 

Preferably, the step of whitening the first electric signal and the 
second electric signal is performed with independent component analysis. 

Preferably, the signal processing method further comprises a step of 
obtaining a signal-to-noise ratio of the observed pulse wave data by performing 
a frequency analysis with respect to the signal component and the noise 
component at each one of the frequency ranges. 

According to the invention, there is also provided a pulse photometer 
In which each of the signal processing methods is executed. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The above objects and advantages of the present invention will 
become more apparent by describing in detail prefenred exemplary 
5 embodiments thereof with reference to the accompanying drawings, wherein: 
Fig. 1 is a block diagram showing the configuration of a pulse 
oximeter according to one embodiment of the invention; 

. Fig. 2A is a waveform obtained through logarithmeric processing ; . 
Fig. 2B is a graph plotted with the amplitude of infi^red light being 
10 taken as a horizontal axis and the amplitude of red light being taken as a 
vertical axis; 

Fig. 3A shows observed (measured) waveforms on v\^ich noise is 
superimposed; 

Fig. 3B is a graph plotted with the amplitude of the infrared light being 
1 5 taken as a horizontal axis and the amplitude of the red light being taken as the 
vertical a^ds; 

Fig. 3C Is a graph obtained by subjecting the graph shown in Fig. 3B 
to processing performed by the transformation matrix A; 

Fig. 3D is a waveform whose noise has been separated or reduced 
20 through the processing; 

Fig. 4 is a view showing the relationship between a signal component 
and a noise component; 

Fig. 5 is a view showing an infrared light spectrum and a red light 
spectrum, t)oth obtained through frequency analysis performed with respect to 
25 the waveforms shown in Fig. 3A; 
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Fig. 6 is a view showing a sp ctrum obtained through frequency 
analysis of the signal compon nt and the noise component shown in Fig. 3D; 

Fig, 7 shows graphs plotted in the same way as that In Fig. 2B for 
every divided frequency ranges, in a case where noise is superimposed on an 
observation signal; 

Fig. 8 shows graphs prepared by the principle shown in Fig. 2B for 
every divided frequency ranges, in a case v\^ere noise is not superimposed on 
an observation signal;. 

Fig. 9 shows a result obtained by determining a gradient for each of 
the frequency ranges shown In Fig. 7, through Independent component 
analysis; 

Fig. 10 shows a result obtained by determining a gradient for each of 
the frequenqr ranges shown in Fig. 8, through the independent component 
analysis; and 

Fig. 1 1 is a flowchart showing a proces^'ng flow performed in the 
pulse oximeter shown in Fig. 1 . 

DETAILED DESCRIPTION O F THE INVENTION 

(How to Separate Signal Component from Noise Component by Utilization of 
Known Light Absorbance Ratio) 

First, the light received by the light receiver 3 provided in the probe 1 
having such a configuration as that shown in Fig. 1 is converted into a voltage 
by the input section 5. 

A component reflecting an optical diaracteristic of pulsation arising In 
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an area subjected to the measurement appears in the received light signal as 
an AC component. 

The signal output from the light receiver 3 is separated and 
demodulated into a waveform of infrared light (IR) and a waveform of red light 
(R) by the demodulator 6, and the waveforms are digitized by the A/D 
converter 7. 

A data sequence of the Infrared light observation signal IR, a data 
sequence of the red light observation signal R,: pulse wave components .of the 
observation signal $, and noise components of the observation signal N are 
respectively taken as: 



As shown in Fig. 4, provided that the gradient of the pulse wave component is 
taken as 4> and the gradient of the noise component is taken as (p, the resultant 
observation signals (IRi, Ri) are expressed as the sum of vectors, Le., a signal 
component Si and a noise component Ni, as follows: 



IR = {IRi; i=K), 1, 2, ...} 
R = {Ri; i=0, 1.2, ...} 
S = {Si; 1=0.1.2. ...} 
N = {Ni; i=0.1,2, ...} 



(13) 
(14) 



(12) 



(11) 



(15) 



Here, defining a matrix W as: 
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w 



f COS( 



(16) 



If a transformation matrix A which is an inverse matrix of W can be 
caused to act on observation data, the signal component S can be separated 
from the noise component N as represented below. 



(17) 



This means that the signal component S is converted Into a horizontal axis and 
the noise component N is converted into a vertical axis. The transformation 
matrix A to be used for implementing whitening is determined as: 



cos^ cos^ 
sin ^ sm<pj 
1 f^sin^ -cos^ 



sin(^-(*)\^-sin$i^ cos^ j 



(18) 



Therefore, we have: 



(19) 



The light absorbance ratio O designates the gradient of a signal 
component The relationship existing between the gradient ^ and the light 
absorbance ratio <p is expressed as O = tan<|>. 

If a period of time is sufficiently shorter than the duration of a pulse 
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rate, an immediately preceding value (e.g., a value obtained 10 msec, 
previously) can l:>e taken as a l^nown value of Alternatively, the light 
ab$ort>ance ratio which has been determined that the noise is not mixed 
therein is used as the known value of O. Still alternatively, a candidate for ^ 
obtained by the technique described later can also be employed. 
Next, 6 is expressed as: 

.e=:<p-4> (20) 

Since <p is unknown, Bis determined as an angle at whidi a norm in the 
direction of the horizontal axis becomes minimum between ^> to (vJ2 - 

When the observation sisr^l is processed through use of the 
transformation matrix A using the thus-obtained 4>, 8, a signal component 
appears along the horizontal axis, and the noise component appears along the 
vertical axis. 

IWore specifically, the determinant expressed as Equation 19 means a 
coordinate transformation from a first coordinate system in which the observed 
signals IR and R are selected as the coordinate axes shown in Fig. 3B, to a 
second coordinate system in which the signal component and the noise 
component are selected as the coordinate axes shown in Fig. 3C. 

f ig. 4 shows observed values (IRi, Ri) at a certain time point t = ti. 
The vector Xi representing the observed values is a synthesis vector of a 
signal component vector Si and a noise component vector Ni. Upon 
performing the above coordinate transformation, the signal component vector 
Si is first operated to be made coincident with the horizontal axis of the second 
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coordinate system. Namely, the vector Si shown in Fig. 4 is rotated clockwise 
by an angle of ^. The value of ^ is obtained through the above relationship 
using the known value of 

The coefficient 6 in Equation 19 corresponds to an operation to make 
the noise component vector Ni coincident with the vertical axis of the second 
coordinate system. Since the waveform shown in Fig. 3B Is a set of the 
observed values, if the values of ^ are made constant for the observed values 
at any of the time points, the values of 8. accordingly vary corresponding to the 
respective time points. In order to uniquely determine the value of 6 as the 
coefficient, disp^sion of values (norm) in the horizontal axis (signal 
component) direction is checked while varying the value of e within a range 
from -4> to (n/2 - 4>). The value of 8 is finally determined as a value making the 
norm minimum. 

Fig. 3C shows a graph obtained by subjecring the graph shown in Fig. 
3B to processing performed by the transformation matrix A. Further. Fig, 3D . 
shows a waveform from which noise has been separated through the 
processing. 

Fig. 5 shows an infrared-ray spectrum (IR) and a red light spectrum 
(R), both being obtained through frequency suialysis performed before 
processing, and Fig. 6 shows a spectaim obtained through frequency analysis 
of a signal component and a noise component performed after corresponding 
processing. 

(How to Calculate Pulse Rate) 

A peak-to-peak interval Ts (Fig. 3D) is obtained from the signal 
waveform extracted through the foregoing conversion processing, whereby a 
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pulse rate = 1/Ts x 60 [beats/min.] can be readily determined. 

Alt matively, the fundamental frequency fs of a pulse wave in relation 
to the spectrum (Fig. 6} obtained from the processed signal through frequency 
analysis becomes dearer than the spectrum (Fig. 5) obtained through 
frequency analysis performed before processing. Hence, a pulse rate can be 
computed by a pulse rate = fs x 60 [beats/min.]. 
(How to Determine Whether Noise Is Mixed into Observation Signal) 

An observation signal is separated into signals each associated with 
one of predetermined frequency ranges whic^ are divided. The division of the 
observed frequency band may be performed for consecutive frequency ranges 
or discrete ranges. 

In each frequency range, an infrared light component is talcen as a 
horizontal axis (or a vertical axis), and a red light component is taken as the 
axis orthogonal to the horizontal axis (or the vertical asds). 

Fig. 7 shows graphs plotted with the infrared li^t component being 
taken as the horizontal axis and the red light component being taken as the 
vertical axis. 

A gradient corresponding to the prindpai component of the graph is 
determined over the entire pulse wave frequency range ("all bands" shown in 
Figs. 7 and 8) and for each frequency range obtained after separation. The 
observation signal is determined to be a single signal or the sum of a plurality 
of signals, on the basis of whether a match exists between the gradients. 

Specifically, a method for determining whether an observation signal 
is a single signal is to determine whether the gradient of the entire pulse wave 
frequency band and the gradients of all the separated frequency range fell 
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within a predetermined range, or whether the gradient of a frequency range 
having a relatively large power and the gradient of the entire pulse wave 
frequency band fall within the predetermined range. 

Specifically, Fig$. 7 and 8 show that the entire pulse wave fr^u^cy 



The gradient of the principal component is described at the top of each graph. 

Fig. 7 shows different values in ranges having a relatively large power, 
such as 0.5 to 1 Hz: 0.54, 3.5 to .4.0 Hz: 0.84, and 4.0 to 4.5 Hz: 0.84. 
Summation of a plurality of signals can be determined to be observed. 

Fig. 8 shows a value of 0.50 is exhibited in a range of 1 to 1.5 Hz, a 
range of 2 to 2.5 Hz, and a range of 2.5 to 3 Hz, each having a relatively large 
power. Hence, the observation signal is determined to be a single signal. 
(Computation of Candidate for Oxygen Saturation) 

The pulse wave signal firequency band is divided into a plurality of 
ranges, as in the case of determination of mixing of noise, independent 
component analysis is performed for each of the thus-separated frequency 
ranges. 

Under the assumption that a longitudinal vector of the observation 
signal is taken as X and a transposed matrix is tal<en as T, the following 
equation is derived from Equations 1 1 and 12. 



band and the range from 0.5 Hz to 6 Hz are divided in increments of 0.5 Hz. 



(21) 




(22) 
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A variance-covariance matrix H of the longitudinal vector X of the observation 
signal is expressed as: 



H^X^^X (23) 

A diagonal matrix having eigenvalues as diagonal elements is taken as A, and 
a matrix formed from eigenvectors is taken as F. Provided that the 
transformation matrix is taken as U, the independent component Y is obtained 
as: 

Y^UX^T'h^f'X (24) 



The independent component Y is dissolved into independent components. 

If signal components are present in the thus-separated frequency 
ranges, the signal components can be converted into a signal 
component-noise component plane such as that shown in Fig. 3C. 

If the matrix formed from eigenvectors is converted through use of an 
inverse matrix U"^ of the transformation matrix U, a vector Is projected onto the 
IR-R plane. The gradient of the vector expressing the signal component is 
obtained as Os, and the gradient of the vector representing the noise 
component is obtained as On. 

Fig. 10 shows results of determination of the gradient of the 
observation signal shown in Fig. 2 on each of the frequency ranges shown in 
Fig. a The light absorbance ratios 1 are understood to be close to each other 
over the respective frequencies (i.e., 0,5±0.02). Fig. 9 shows results obtained 
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by processing the graphs shown In Fig. 7, and a result showing *s is obtained 
within a range including fs shown in Fig. 6, The oxygen saturation Sp02 is 
obtained as a function of <(>s; that is, Sp02 = f (Os). 

Operation of the processor of the device block shown in Fig. 1 will 
5 now be described in detail with reference to Fig. 1 1 . 

When measurement is commenced, digital signals obtained as a 
result of observation of infrared light and red light are subjected to preliminary 
processing (step S1}. Here, the preliminary processing refers to logarithmic 
computation or pulse wave band filtering on the basis of the Lambert*-Beer law. 
10 The signals are divided into predetermined frequency ranges (e.g.. In 

increments of 0.5 Hz) through use of band-pass filths, thereby determining a 
principal component of light absorbance ratio (i.e., the gradient of the principal 
component on an IR-R plane) for each frequency range (step S2). 

\AAiether a match exists between the gradients det^mined for the 
15 respective frequency ranges is determined, thereby determining a candidate 
for a light absorbance ratio of a pulse wave (step S3). A plurality of 
candidates may be acceptable. 

A determination is made as to whether the candidate (signal source) 
is single or plural (step 84), on the basts of the results obtained in the step S2. 
20 In a case where the determination is single, oxygen saturation is 

determined from the principal component light absorbance ratio, and the 
fundamental cyde is determined from the result obtained in the step S1. 
Thus, a pulse rate is computed (step S5). 

In a case where the determination is plural, at the step S6, affine 
25 transformation is performed through use of the candidates for the light 
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absorbance ratio obtained in the step S3 or the result of the light absorbance 
ratios corresponding to an immediatety-preceding Sp02 valu (sufficiently 
shorter than the pulse rate: for example, 10 msec faster), or a light absorbance 
ratio which has been determined that the noise is not mixed therein. Equation 
19 is used for the transformation matrix, and the matrix Is rotated with 
reference to ^ (i.e., a value corresponding to a light absorbance ratio), and 9 is 
determined as an angle at which a norm value in the direction of the ^ axis 
becomes minimum. .... 

In order to idmtrfy a pulse rate, the fundamental cycle Ts (see Fig. 
3D) is determined from a cycle determined along the time axis of the signal 
obtained in the direction of 4> determined in the step S6, or the frequency 
analysis is performed to identify the fundamental frequency fs (see Fig. 6). A 
signai-to-noise (S/N) ratio Is determined for each frequency, and the 
thus-determined ratio is taken as one index for evaluating the reliability of the 
obtained light absorbance ratio (step S7). 

On the other hand, time domain data for the respective frequency 
ranges, which are divided at the step S2, are used to perform independent 
component analysis (step S8). 

The gradients at the IR-R plane of the eigenvectors obtained in the 
step S8 are determined, and it is checked whether a match exists between the 
frequency ranges, to thus determine a candidate for a light absorbance ratio 
(step S9). A plurality of candidates may be acceptable. 

Outputs obtained In the steps S3, S5, S7, and S9 are compared with 
the light absorbance ratio and bands obtained in the past. Current oxygen 
saturation is determined from a common item, consistency, and reliability 
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(output from the step S7). The oxygen saturation is averaged, as required, 
and the thusnaveraged oxygen saturation is output (step S10). 

Wlien measurement is performed continuously, processing returns to 
the step S1, and processing is iterated. If measurement Is not p^formed 
continuously, measurement is terminated. 

According to the invention, the pulse wave signal is separated from 
the noise signal through use of a known light absorbance ratio, thereby 
acquiring a pulse wave signal whose noise has been.reduced. 

affine transformation Is used for utilization of a known light 
absorbance ratio, thereby whitening an observation signal. The pulse wave 
signal is separated from tlie noise signal. Thus, a pulse wave signal whose 
noise has been reduced can b& acquired. 

The pulse wave signal frequency band is divided into a plurality of 
frequency ranges. The prindpal component of the light absorbance ratio is 
determined for each frequency band. By Judging whether a match exists 
between the light ab$ort>ance ratios, a determination can be made as to 
virfiether noise is mixed into the observation signal. 

The afRne transformation utilizing a known light absorbance ratio 
enables determination of a first rotation utilizing a known light absorbance ratio 
and a coefficient suitable for reducing noise. Thus, a pulse wave signal 
whose noise has been reduced can be acquired. 

A pulse can be detected from the pulse wave signal whose noise has 
been reduced, or a pulse rate can be determined through frequency analysis of 
the pulse wave signal. 

Dividing the pulse wave signal frequency band into a plurality of 
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frequency ranges enables a reduction In the number of signal sources. 
Further, performing whitening operation on a per-frequency-range basis 
enables deteanination of a light absorbance ratio of the pulse wave signal and 
a light absorbance ratio of the noise signal. 

Since the pulse wave signal is separated from the noise signal, an 
S/N ratio is obtained for each divided frequency range, and the thus-obtained 
S/N ratio can be utilized as an index to be used for evaluating a signal. 
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